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C~| . Abstract 

Nonlinear dynamics of a bouncing ball moving in gravitational field 
and colliding with a moving limiter is considered. Several simple models 
of table motion are studied and compared. Dependence of displacement of 
the table on time, approximating sinusoidal motion and making analytical 
^ ' computations possible, is assumed as quadratic and cubic functions of 

, time, respectively. 
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1 Introduction 

\D _ 

| Vibro-impacting systems belong to a very interesting and important class of 

nonsmooth and nonlinear dynamical systems [TJ [3J [31 S] with important tech- 
nological applications [5J [HI Ul IE] • Dynamics of such systems can be extremely 
complicated due to velocity discontinuity arising upon impacts. A very charac- 
teristic feature of such systems is the presence of nonstandard bifurcations such 
as border-collisions and grazing impacts which often lead to complex chaotic 
motions. 

The Poincare map, describing evolution from an impact to the next impact, 
is a natural tool to study vibro-impacting systems. The main difficulty with 
investigating impacting systems is in finding instant of the next impact what 
typically involves solving a nonlinear equation. However, the problem can be 
simplified in the case of a bouncing ball dynamics assuming a special motion 
of the limiter. Bouncing ball models have been extensively studied, see [9] and 
references therein. As a motivation that inspired this work, we mention study 
of physics and transport of granular matter [B] . A similar model has been also 
used to describe the motion of railway bogies [7] . 

Recently, we have considered several models of motion of a material point 
in a gravitational field colliding with a limiter moving with piecewise constant 
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velocity [TOJ HU HH EH H] • In the present paper more realistic yet still simple 
models are considered. The purpose of this work is to approximate sinusoidal 
motion as exactly as possible but still preserving possibility of analytical com- 
putations. 

The paper is organized as follows. In Section 2 a one dimensional dynamics 
of a ball moving in a gravitational field and colliding with a table is considered 
and Poincare map is constructed. Dependence of displacement of the table 
on time is assumed as quadratic and cubic functions of time, respectively. In 
Section 3 bifurcation diagrams for such models of table motion are computed 
and compared with the case of sinusoidal motion. 

2 Bouncing ball: a simple motion of the table 

We consider a motion of a small ball moving vertically in a gravitational field 
and colliding with a moving table, representing unilateral constraints. The ball 
is treated as a material point while the limiter's mass is assumed so large that 
its motion is not affected at impacts. A motion of the ball between impacts is 
described by the Newton's law of motion: 

mx = —mg, (1) 

where x — dx/dt and motion of the limiter is: 

y = y(t), (2) 

with a known function y. We shall also assume that y is a continuous function 
of time. Impacts are modeled as follows: 

%(n) = y(n), (3) 
±{rt)-y(n) = -R(xfc)-y(T$), (4) 

where duration of an impact is neglected with respect to time of motion between 
impacts. In Eqs. <j3j> , (jlj n stands for time of the i-th impact while x~ , xfaie 
left-sided and right-sided limits of ii (t) for t — > n , respectively, and R is the 
coefficient of restitution, < R < 1 [5]. 

Solving Eq. (fTJ) and applying impact conditions ©, (@| we derive the 
Poincare map }15j : 

7 y(T m ) = jY (Ti) - A? +1 + A i+x Vi, (5a) 
V l+1 = -RV i + 2RA i+1 + 1 (l + R)Y{T i+1 ), (5b) 

where A^ + i = Ti + \ — Ti. The limiter's motion has been typically assumed 
in form Y S (T) = sin(T), cf. [11 and references therein. This choice leads to 
serious difficulties in solving the first of Eqs. ([5]) for Tj+i, thus making analytical 
investigations of dynamics hardly possible. Accordingly, we have decided to 
simplify the limiter's periodic motion to make (|5a[) solvable. 
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In our previous papers we have assumed displacement of the table as the 
following periodic function of time: 



Y(T) 



if, < f < h 

h<f<l 



l-h 1 - 



1 



(6) 



with T = T — \ T\ and < h < 1, where [x\ is the floor function - the largest 
integer less than or equal to x. In this work the function Y(T) is assumed as 
quadratic Y q and two cubic functions of time, Y Cl and Y C2 . More exactly, these 
functions read: 



Y C2 (T) 



Y q (T) 
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< T < 
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Y C1 (T) = 12V3T (f - i 



T - 1 



< T < 1 



< T < 1 



(7) 
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(32?r - 128) f 3 + (-16tt + 48) f 2 + 2ttT 
(128 - 32tt) f 3 + (-144 + 32tt) f 2 + (48 - IOtt) f - 4 - 
(128 - 32tt) f 3 + (-240 + 64tt) f 2 + (144 - 42tt) f - 28 
(32tt - 128) f 3 + (336 - 80tt) f 2 + (-288 + 66tt) f + 80 - 
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- 9tt 
18tt 



< T < \ 

1 < T < i 
4 — J ^ 2 

1 < J 1 < 3 

2 — J ^ 4 



< T < 1 



see Figs. 1 ,2, 3. 



(9) 



Y(T) 




Figure 1: Displacement of the table, Y q {T) (black) and Y S (T) (red). 

The function Y q , consisting of two parabolas, and its first derivative are 
continuous, however its second derivative is discontinuous at T = 0, ^, 1. 
The function Y Cl is smooth but provides a poorer approximation to Y s (T) = 
sin(27rT) then Y q . We have included this function because it is smooth and is 
the lowest-order polynomial approximating Y s on the unit interval [0, 1]. 
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Figure 2: Displacement of the table, Y Cl (T) (black) and Y S (T) (red). 



The smooth function Y C2 consists of four cubic functions and provides the 
best approximation to Y s . Let us note that for all these functions, Y q , Y Cl , Y C2 , 



equation for 




Figure 3: Displacement of the table, Y C2 (T) (black) and Y S (T) (red). 

Our model consists thus of equation (J5j> and one of Eqs.Q, (jHJl, © with 
control parameters R, 7. Since the period of motion of the limiter is equal to 
one, the map ([5]) is invariant under the translation Tj — > Tj + 1. Accordingly, 
all impact times Tj can be reduced to the unit interval [0, 1]. 
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3 Comparison of bifurcation diagrams 



We have computed bifurcation diagrams to study dependence of dynamics on the 
model of motion of the table. Let us recall that in the case of displacement of the 
table described by piecewise linear function Y defined in Eq.® the bifurcation 
diagram differs significantly from that computed for sinusoidal displacement Y s . 




Figure 4: Bifurcation diagram, Y(T) = Y q (T) 




Figure 5: Bifurcation diagram, Y(T) — Y Cl (T) 



In the h — > 1 case only manifolds of periodic solutions were found to exist 
and there is no chaotic dynamics [TT] while for h ^ 1, classical attractors exist 
as well, but only one period doubling on route to chaos via corner bifurcation 
was reported [13] in contradistinction to the case of sinusoidal displacement of 



5 



the table where full period doubling scenario is generic [15] . 

The bifurcation diagram for the displacement function Y q is shown in Fig. 
4. There are one 2-cycle (ending after one period doubling) and six fixed points 
(first three are shown) which become unstable at some 7 cr where chaotic bands 
appear. There are no manifolds of periodic solutions. The bifurcation diagrams 
for functions Y C1 and Y C2 are presented in Figs. 5 and 6, respectively. It can be 
seen that full cascades of period doubling with transition to chaos are present 
in Fig. 5 but it seems that some bifurcation paths end abruptly. Moreover, 
periodic states loose stability in different order than for Y s . 




Figure 6: Bifurcation diagram, Y(T) — Y C2 (T) 




Figure 7: Bifurcation diagram, Y(T) = sin (2irT) 
We note finally that the bifurcation diagram computed for the displacement 



function Y C2 , cf. Fig. 6, is very similar to the bifurcation diagram for the sine 
function Y s , see Fig. 7. 

4 Summary and discussion 

We have constructed several simple models of table motion in bouncing ball 
dynamics in order to approximate sinusoidal motion as exactly as possible. In 
conclusion we can state that dynamics of the model based on Eqs. (J3J) , © cor- 
responds well to dynamics with table displacement given by Y s (T) = sin (2irT). 
Moreover, equation for time of the next impact (|5aj) is a third-order algebraic 
equation in Tj+i for Y (T) = Y C2 (T) and thus analytical computations are pos- 
sible. We are going to investigate all these models in our future work. 
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